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Dynamical response of nanomechanical cantilever structures immersed in a viscous fluid is impor- 
tant to in vitro single- molecule force spectroscopy, biomolecular recognition of disease-speciflc pro- 
teins, and the detection of microscopic dynamics of proteins. Here we study the stochastic response 
of biofunctionalized nanomechanical cantilevers beam in a viscous fluid. Using the fluctuation- 
dissipation theorem we derive an exact expression for the spectral density of the displacement and 
a linear approximation for the resonance frequency shift. We find that in a viscous solution the 
frequency shift of the nanoscale cantilever is determined by surface stress generated by biomolecular 
interaction with negligible contributions from mass loading. 

PACS numbers: 81.07.-b,45.10.-b,83.10.Mj,82.70.Dd 



From single molecule force spectroscopy 1] to biomolec- 
ular recognition of disease-specific proteins such as can- 
cer antigens 0], micron-sized cantilevers have proved to 
be fundamental to the ultrasensitive detection of small 
forces. Usually, forces are detected by measuring the de- 
flection of the cantilever. In the dynamic case, the shift 
in the resonance frequency of the cantilever is used to 
infer the magnitude of the force. Micromachining tech- 
niques now enable commercial production of such can- 
tilevers with dimensions on the order of 100 ^m as well 
as their routine use in force spectroscopy. 

Decreasing the cantilever dimensions to sub-micron or 
nanometer scales increases the resonance frequency to the 
megahertz-gigahertz range. The resultant increase in the 
dynamic range and the measurement speed can provide 
a better tool for probing single molecules. This could be 
also used for more sensitive bioimaging techniques and 
monitoring real-time binding kinetics of ligand-protein 
binding as well as the energy landscape of the molec- 
ular bonds at their true characteristic time scales. For 
biomolecular recognition in a viscous fluid, force sensitiv- 
ity can be increased by decreasing the effective viscous 
damping. Nanoscale cantilevers are hence expected to 
have dramatically enhanced force sensitivity as smaller 
cantilevers have lower viscous damping. 

In spite of the obvious importance of nanomechan- 
ical cantilevers for ultrasensitive in vitro force detec- 
tion, there is no widely-accepted description that relates 
resonance frequency change to concentration or mass- 
loading, over the entire range of viscosity, relevant to 
biomolecular recognition in viscous fluids. In this Letter, 
we derive an exact expression that takes into account the 
hydrodynamics of the beam in a continuum approxima- 
tion with axial loading. We find that, in air, the first 
mode of the nanoscale cantilever is more sensitive to sur- 
face stress than the higher order modes. The frequency 
shift is primarily determined by surface stress in the first 



mode, where the mass-loading effects become relevant for 
higher order modes. More importantly, we find that, in 
a viscous solution such as water, frequency shift is domi- 
nated by surface stress and not mass loading, as generally 
expected. 

Model: In standard fiuid dynamics, an inviscid model 
is typically a valid assumption. However, nanoscale can- 
tilevers have Reynolds numbers less than one, therefore 
viscous effects become dominant 0|. Proper inclusion of 
viscous effects is characterized as being in the Stokes- 
Purcell Regime of fiuid fiow. 

Neglecting rotatory inertia, shear deformation and in- 
ternal damping, the equation of motion for the defiection 
y{x,t) of a beam with length L, width h and thickness 
d, immersed in a fluid at temperature T and loaded by a 
constant axial force S, is given by 

E and / are the Young's modulus and moment of inertia 
of the (coated) beam, respectively. fh{x,t) is the hydro- 
dynamic loading due to the motion of the fluid around 
the beam and fthix,t) is a Langevin-type force per unit 
length responsible for the thermalization of the beam. 
The linear mass (mass per unit length) of the system n{x) 
consists of the linear mass of the beam /if,, and the linear 
mass of the trapped biomolecules fJ-i{x). The axial load 
S introduced in describes the mutual interaction of 
biomolecules adsorbed on the beam The boundary 
conditions for Eq. (Q) are given by 2/(0, t) = y'(0,t) = 0, 
y"{L,t) = and Ely'" (LA) = Sy'{L,t), where primes 
denote spatial derivatives |l2l |. 

At higher concentrations, the biomolecules form a uni- 
form layer with mass fii and thickness h so that 

fi{x) — fib + fJ-i ^ constant. (2) 

The mutual interaction of biomolecules within the layer is 
modeled by taking into account the stress a they generate 
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on the coated surface of the beam. As shown in ref.[3, 
this stress enables the bending of the sihcon-nitride mi- 
crocantilevers with length to thickness ratio L/d, ranging 
from 10^ to 10^. The resulting static deflection, on the 
order of a few tenths of microns, is related to surface 
stress by Stoney's formula 10]. However, for the sili- 
con nanomechanical cantilevers under investigation here 
{L/d ~ 50), Stoney's formula typically yields angstrom- 
level bendings (lO^^L). So these nanocantilevers remain 
almost straight under the influence of surface stress. Nev- 
ertheless, as shown in this stress induces an effective 
axial load, S — aL, that must be included in Eg . ffl . In 
vacuum, such a model has been studied in ref.jjj]- 

In addition, biomolecular interaction on the surface re- 
sults in an effective Young's modulus of the layer Ei [see 
ref.0]: 

EI = Ebh + Eih ~ Ebbd^/12 + Eihhd^/A. (3) 

Ebh and Eili are the respective bending rigidities of the 
beam and the layer (the last equality holds when /i <C d) . 
In the opposite limit, biomolecules with mass m sparsely 
scattered over the beam at locations Xi, result in 

n{x) ^ ^lb+^Jil{x) ^ ^ib+m'^5{x-Xi), Xi^[Q,L\. (4) 

If the average spacing is large compared to their size, 
their mutual interaction is negligible and 5 = 0. Consid- 
ering that their presence does not substantially affect the 
moment lb of the beam, the bending rigidity of the whole 
system is the same as for an unloaded beam, EI — Eblb- 
Equations of Motion: To solve equation , we expand 
the deflection y{x,t) and the force densities fh{x,t) and 
/th(a;, t) in terms of the modes of the bare beam, defined 
as the beam without the added mass {fii {x) = 0) though 
it includes the tension S = aL: 

oo 

y(a;,i) = ^2/„(t)0„(C), e = x/L. (5) 

Similar expressions hold for the force densities. The 
eigenmodes 4>niO satisfy the following conditions: 

0„(O) = = = 0, = a-^Ul)- (6) 

The self-adjointness of © makes the modes orthonormal: 

Jo4'n{04'i{(,) = Sn,i. From Eq. the eigenvalues /3„ 
are the successive positive roots of 

1 + (1 + £^j) cosh A+ cos + e„ sinh A+ sin A^ = (7) 

where e„ = a/(2/3^) and A^ — /3n[\/l + ± e„]5. Note 
that, for e„ — a/ (2/3^) <C 1, Eq. Q reduces to the usual 
clamped-free equation, 1-1- cos /3„ cosh/3„ = 0. As /3„ oc n 
for large n, eigenvalues for which n ^ ^/a are essentially 
independent of the surface stress. 



Let ^{x) = ^J,b + Ih{x)- Using Eq. Q reduces to 

oo 

Myn{t)^Kyn{t)+Y,'^^oyAt) = FnAi)+Fn,th{t). (8) 

M = pLbL is the mass of the beam and J^n,h(th)(^) — 
L/„ The effective stiffness of mode n is 

K - EIl3fjL^ (9) 

The real and symmetric matrix $ has components 

= l/ M,(e'i)0„(r)0,(r)de'. (10) 

If fJ.i{x) = fii, then €> = M;l, where 1 is the identity 
matrix and Mi = L^i is the layer mass. Eqs. (O decouple 
and the mass of modes ?/„ becomes the total mass of the 
system, M + Mi. But non-uniform mass distributions as 
in Eq. couple the bare modes of the beam. 

Taking the Fourier transform of and using the ex- 
pression for the hydrodynamic force 

F„,h(w) = Mfiu^T{io)ynioj), (11) 

where Mf — jLpfb^ is the mass of the fluid loading the 
beam, and T{uj) = Tr{Lo) + iTi{uj) is a complex "hydro- 
dynamic function" discussed in detail in we obtain 

A(u;)|y(c.)) = \F,,,{lu)). (12) 

Kets \v) are column vectors with components Vi, i e N. 
The nonhermitian matrix A (a;) is given by 

A(a;) = Ao(a;)-w2* (13) 

where 

Ao(tj),y = [fc„ - Lu'^iMfTiuj) + Mn)]5nj 
^nj = *nj(l - Snj) ; M^ = M + $„„. 

Spectral Densities: As the dissipative (complex) part 

of the hydrodynamic function is frequency-dependent, we 

apply the generalized fluctuation-dissipation theorem 

to derive the power spectrum matrix of the stochastic 

— s 

forces F„^th, Sp{u}) = |Fth(w))(Fth(w)| (the over-line 

denotes thermal averaging, the superscript s refers to the 

spectral density and (i^thl is the hermitian conjugate of 

lAh)): 

kT 

Sp{Lo) = — (A^fiuj) - A(a;)) = 2kTMfLur,{uj)A (14) 

where k is the Boltzmann constant and T the tempera- 
ture. In components, this yields 

Fn,th(w)-F;*th(w') = 2kTMfUjT,{uj)SnpS{uj - uj'). (15) 
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Notice that this expression does not depend on It 
is the same as for a bare beam. As seen above, the 
stochastic forces acting on distinct modes are uncorre- 
lated. Nevertheless, their power spectrum cannot be as- 
sumed to be constant, contrary to the assumption in ref. 
Q , as the dissipative part of the hydrodynamic function 
is frequency-dependent. Eq. (|15() is the generahzation 
of the expression derived by Paul and Cross for a single 
cantilever mode 0- Now, inverting l|12|l . we obtain 

=x('^)|AhM), x(a;) - A-i(a;) (16) 

and \yiLo)){yiu;)\ = x(u;)|Fth(a;))(Fth(w)|x^(u;). Spec- 
tral averaging the latter and using H14|) , we find the power 
spectrum matrix of the deflection modes 



(17) 



Introducing with components (f>niOj the Fourier 

transform of the deflection (0 reads y{x,uj) — {ip^\y{Lu)) 
and using (|17|l . we find its spectral density to be 



\y{x,uj)\^ ^2kTMfCur,{iu){cl>^\xiuj)xH^)\k)- (18) 

The total mass of the particles trapped on the beam Mi is 
small compared to the mass of the beam. This justifies 
treating a>^# perturbatively provided its elements stay 
small compared to the diagonal elements of Ao(u; 
is indeed the case provided MfTi{Lu) ^ AMi 
inverting (|13|l . we obtain in first order in $ 



This 
By 



XX^ = Xoxl + (*Xo + xl^) xl + 0(*') (19) 
where Xo = -'^o ^- Reinstating in ifT^ . we finally get 



{x,u})\ ^2kTMfUjT,{uj) 



E 



M2M. 



oo oo 



+2^^ EE 



0(*2) 



(20) 



In this expression, S, = x/L, $„p is given by (|14(l and the 
quantity A„ = Ao(a;)„„/M„ reads 



A. = -c^'(i + A„r(c^)) 



(21) 



where u;„ = yJkn/Mn is the frequency in vacuum, A„ = 
Mf/Mn and 9?(A„) is the real part of An- An expression 
similar to (|20|l can be derived for the slope of the deflec- 
tion provided (f>j{^) is replaced by (f>j{£,) and the overall 
prefactor is divided by L^. 

Expression 1)2011 is valid for any mass distribution /i;(a;) 
along the beam. For a uniform layer with linear mass 
^J'l, ^np — Lfiidnp, and then $„p — 0. All modes have 
the same effective mass, M„ = M + Ljii, and are de- 
coupled. Reinstating in H20|l . the second term vanishes 
and we obtain the exact spectral density of a composite 



8 
7 

6 l- 
5 
4 
3 
2 
1 



Exact 
Sader 



/ = uj/2n (MHz) 

FIG. 1: Spectral density of the deflection slope for a 
silicon nanobeam in air and water (solid lines) compared to 
Sader's expression Q (dashed lines). 



beam consisting of the original beam plus the layer. For 
molecules trapped on the beam at positions ^i, the mass 
profile given in Eq.Q leads to <i>„p ^ mj^i 0n(CO0p(Ci) 
and Ea. (|20|l is valid up to first order in m provided the 
frequency satisfies Mfri{uj) ^ 4M;. Interestingly, if we 
assume N molecules to be randomly scattered along the 
beam in a uniform way and average $„p accordingly, we 
find (($„p)) ^mY^iJ (l^n{(,i)4>p{ii)d£,i = Nm5„p- As the 
total mass of the trapped molecules is small compared 
to the mass of the beam, in the first approximation, the 
average spectral density is the same as the spectral den- 
sity of their average mass distribution — i.e. the spectral 
density of a uniform layer of mass Nm. 

In Fig. n we compare the bare beam ($ = 0) spectral 
density of the deflection slope at the tip of a rectangular 
silicon nanocantilever {E = 160 Gpa, p = 2.33 x 10'^ 
kg/m^) to Sader's result 4] in air and water. The 
beam dimensions are d x b x L = 0.2 x 0.2 x 10(/im). 
From at T = 27° C, the viscosities are r/air = 
1.86 X 10-^ Tywatcr = 8.59 x 10"'' (kg/m/s) and the den- 
sities, Pair = 1-18, Pwater = 997 (kg/m^). Afihough dif- 
ferent, Sader's formula can be shown to reduce to H20|l 
provided A|r(cj)| <C 1. This explains why the results are 
very similar in air (A|r(a;)| ~ 0.06 at resonance) while 
they start to differ in water (A|r(tt')| ~ 25 at resonance). 

Frequency Shift: As stated earlier, when trapped 
molecules form a uniform layer, the exoct spectral density 
of the beam deflection is given by 



(X,UJ 



2kTMfUjT^{u 



OO 

E 



|A„P' 



(22) 



where Mi is the mass of the layer and where An is given in 
Eq.(ini) with M„ = M + Mi. When the peaks of Ea. lf^ 
are sharp enough, the hydrodynamic function is almost 
constant in their vicinity and the resonant frequency sat- 
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FIG. 2: Relative frequency shift (in %) in air and water vs the 
dimensionless stress, a = aL^ /EI . Mi = A//100 and Ei — 0. 
Other parameters are the same as in Fig. Sohd hne: exact 
result from Eq. 1221 . Dotted line: linear approximation from 
Eqs. and (1^. 



isfi.es the self-consistent equation 



IR 
3~ 



^4R^ + 3/2 



(23) 



where R = AIn + MfTr, I = MfTi and Tr^i = ^r,ii'-^R,n)- 
From the expression (|23|l . the mass and stiffness varia- 
tions due to the layer, SM and (5fc„, induce a relative 
frequency shift between a bare and a loaded beam: 



1 



^R, 



dlnf 



Sin/ e,. 
dM 



(24) 



Here, 6M = Mi and its prefactor in H24|) takes into ac- 
count dissipative [MfTi) and fluid mass loading {MfTr) 
effects. According to (O, the two contributions to the 
stiffness 5fc„ come from the bending rigidity, -Efc/f, — *■ 
Eijlt+Eili , and from the surface stress through the eigen- 
value /3„(a). From JSJ and Q, we find 



Sk^ ^ 3m 2r„^„ + /3o,„(r„ + 1„) 

fc„ dEb Pln{in-T^) 



-a. 



(25) 



where /3o,n is the nth root of cosh(/3) cos(/3) -f 1 = and 
where T„ = tanh/3o,n, = tan/3o.n- The last term of 
(I25|l has been obtained from in perturbation, ft is 
valid when a <^ tt? and vanishes as n oo. 

Using the same data as in Fig. ^ we display in Fig. [21 
the relative frequency shift in air (left) and water (right) 
versus the dimensionless surface stress a. The exact 
shift is evaluated from the spectral density (|22(l and 
compared to its linear approximation (|24|l . The layer 
mass has been arbitrarily fixed to 1% of the beam mass 
and El set to zero, hence the negative offset observed 
in air at a = 0. For typical values of the surface 

in ref. 0), 
l/M. Then, 

dlnf /dM ~ —l/M, d\nf/dujR_n ~ 0, and we recover 



stress, a ^ 10 ^J.m ^ (see Wu et al 
\a\ < 1. In air, X\T{lurji)\ < 1, and / 



the usual frequency shift for a linear oscillator in vac- 
uum. As seen on the left panel, the first peak is the most 
sensitive to a. The deviation of the data from the lin- 
ear result H25|l indicates that the condition a <C with 
n — 1 becomes violated. This effect disappears for the 
second and third peaks that are less sensitive to a. In 
water (right panel), a single broad peak occurs. Eq. I|23() 
loses its accuracy but the frequency shift (|24(l derived 
from it is still acceptable. The contribution of d\nf/dM 
becomes negligible while dlnf/dujR^n ~ —dlnTi/dujR^n 
becomes important, hence the increase in the slope of 
the relative frequency shift versus a in water compared 
to air. 



In conclusion, we treat the stochastic response of bio- 
functionalized nanomechanical cantilevers with a gener- 
alized fluctuation-dissipation relation. In a viscous fluid 
like water, the resonance frequency shift for a continuous 
distribution of biomolecules on the cantilever surface is 
dominated by surface stress rather than mass loading. 
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